A unitary correlation operator method 



H. Feldmeieri , T. Neff^ , R. Roth^ and J. Schnack^ 

Gesellschaft fiir Schwerionenforschung mbH, 
Postfach 110 552, D- 64220 Darmstadt & 
Technische Hochschule Darmstadt 



Abstract 

The short range repulsion between nucleons is treated by a unitary correlation op- 
erator which shifts the nucleons away from each other whenever their uncorrelated 
positions are within the replusive core. By formulating the correlation as a transfor- 
mation of the relative distance between particle pairs, general analytic expressions 
for the correlated wave functions and correlated operators are given. The decom- 
position of correlated operators into irreducible n-body operators is discussed. The 
one- and two-body-irreducible parts are worked out explicitly and the contribution 
of three-body correlations is estimated to check convergence. Ground state energies 
of nuclei up to mass number ^4 = 48 are calculated with a spin-isospin-dependent 
potential and single Slater determinants as uncorrelated states. They show that the 
deduced energy- and mass-number-independent correlated two-body Hamiltonian 
reproduces all "exact" many-body calculations surprisingly well. 
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1 Introduction 



In many-body physics the most convenient trial states are products of single- 
particle states which for identical bosons or fermions are symmetrized or anti- 
symmetrized with respect to particle exchange. An adequately selected set of 
those product states spans a subspace of the Fock space called model space. 
In the following we shall use the term product state for bosons and fermions 
as well and proper statistics is implied without mentioning. 

For particles with a rigid internal structure like atoms or nucleons the two- 
body interaction between the particles contains a short-ranged repulsive core. 
For this situation product states are not appropriate because a numerically 
tractable number of them can usually not represent the strong depletion of 
probability in the two-body density distribution at short distances between 
the particles. Or in other words, the true eigenstates of the exact Hamiltonian 
contain components which lie outside the model space in the so-called Q- 
space. 

There are in general two types of methods to remedy the incompatibility of 
the Hamiltonian with the trial states which span the model space [1]. The 
first is to adapt the Hamiltonian by replacing the interaction with an effective 
one which takes the scattering between model space and Q-space into account. 
The Brueckncr G-Matrix formalism is of this type [2] . The method of effective 
operators allows to calculate the whole spectrum of low lying states including 
transitions between them. 

The second way is to introduce the short-range repulsive correlations into the 

many-body state by applying correlation operators [3] onto the trial prod- 
uct state. In the simplest case it is a multiplication with Jastrow correlation 
functions [4]. This approach is mainly used to find ground-state energies by 
minimizing with respect to the parameters contained in the correlated trial 
state. Spectra and transitions cannot be calculated easily in this scheme. 

In this article we propose the Unitary Correlation Operator Method (UCOM), 
which to a certain extend combines the advantages of the above mentioned 
schemes. It correlates the product state, | $ ) , by means of a unitary correlator, 
C. The identity 

(^1 ii" 1^') = ($1 C^i/C 1$') = ($1 C-^i/C 1$') = ($1 ^1$') (1) 

shows that one may either speak of the original "untamed" Hamiltonian H 
sandwiched between correlated states | \t' ) = C | $ ) or of a correlated Hamil- 
tonian H C represented in the model space of the uncorrelated product 
states I $ ) . The task of the correlator C is to "tame" the short-ranged repul- 
sive part of H by means of a unitary transformation. This general idea exists 
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for a long time [5-7] but has not been pursued very much. 

When the correlator C is applied to the Hamiltonian the method can be re- 
garded as a pre-diagonalization such that large matrix elements connecting 
model space and Q-space are removed. H = H C can then for exam- 
ple be diagonalized in a shell model configuration space in order to describe 
long-range many-body correlations in low lying eigenstates. In this sense it is 
similar to the method of effective operators, although there a projection onto 
the model space is used rather than a unitary transformation. 

When C acts on the uncorrelated state | $ ) one has more the picture of 
a trial ansatz for the many-body state as for example used in Jastrow-type 
methods. However, the advantage of a unitary C is that the norms ( \t' | \t' ) = 
( $ I C I $ ) = ( $ I $ ) of correlated and uncorrelated states are identical 
so that there is no need to develop expansion methods when calculating the 
denominator in expectation values like ( $ | C^HC | $ )/( $ | C^C | $ ) [8]. 

We propose to write C as a generalized shift operator, which moves two parti- 
cles away from each other whenever they are within the classically forbidden 
region of the repulsive core. In section 2 we discuss the many-body nature of 
the correlator C and of the correlated operators and give analytic expressions 
for C and the correlated Hamiltonian in terms of a coordinate transformation 
of the relative distance. In the two-body space the resulting correlated Hamil- 
tonian has no repulsive core anymore, instead a non-locality occurs at short 
distances which is written in terms of a quadratic momentum dependence with 
a reduced mass which is a function of the distance. The eigenstates, bound as 
well as scattering states, of the correlated Hamiltonian coincide with the ex- 
act eigenstates of the original Hamiltonian outside the range of the correlator, 
therefore the two interactions are completely phase-shift equivalent [9,10]. 

For the many-body system our aim is to devise a correlator which does not de- 
pend on energy and particle number so that the same correlated Hamiltonian 
H can be used for many nuclei. In section 4 results will be given for mass num- 
bers ranging from A = 2 to A = 48. The main problem consists in separating 
the variational degrees of freedom residing in the correlator C from those in 
the product state | $ ) such that the irreducible three- and higher-body parts 
of H remain small. These aspects are discussed in section 2.4 and in section 
3, where the concept is applied to spin-isospin-independent forces and small 
nuclei first. The convergence of the expansion of the correlated Hamiltonian 
into one-, two- and three-body parts is discussed in terms of a smallness 
parameter which is the product of the density times a correlation volume. 

Then in section 4 the correlator is extended to spin-isospin-dcpcndent forces 
which are more realistic for heavier nuclei. We calculate ground state energies 
up to ^^Ca and compare them to existing "exact" results. The agreement is 
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astonishing considering that the uncorrelated many-body states are single 
Slater determinants. 



2 The concept of a uniteiry correlator 

In this section the advantage of representing improved trial states by a unitary 
transformation of product states will be elucidated in a quite general way. The 
following subsections will explain in detail the explicit form of the correlator C 
and the expansion of a correlated operator into n-body operators, in particular 
the expansion of the correlated Hamiltonian. 

Throughout this section we use the convention that | \E' ) or \ip) stand for 
correlated states and | $ ) or | ) for uncorrelated states. The states | ) 
denote a basis set of many-body product states (symmetrized or antisym- 
metrized if necessary) . 

The eigenstates | n ) of the many-body eigenvalue problem 

H \ n) = En\n) (2) 
can be written as 

\n)^C[H,^ (3) 
where C [H, $] is the unitary operator 

C[i/,$]=E \k){^u\ (4) 

k 

which diagonalizes the Hamiltonian H. 

The idea is to devise a unitary correlation operator C which approximates 
C[i7, $] for low energies, i.e. for states |A;) with Ek in the vicinity of the 
ground state. Prom (4) it is clear that C depends on the choice of the product 
states I ) and on the Hamiltonian H via its eigenstates \k). 

Since there is a certain amount of freedom in choosing the single-particle states 
which build the product states | $„ ) there is no need to be so ambitious to 
really approximate C[H, $] in the sense of an operator equality C — C[H, $]. 
One rather divides the burden among C and | $„ ) in such a way that | $„ ) 
or linear combinations of a small number of | $n ) describe the long-range 
correlations (mean field, deformations) while C takes care of the short-range 
repulsive correlations. 
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In other words the physical task of C is to provide a mapping between the 
Hilbert space spanned by low lying correlated states | k ) and the Hilbert 
space of uncorrelated product states | $fc ) which do not contain (repulsive) 
short-range correlations. 

One main advantage of a unitary correlator C is that the correlated state 

|*) = C|$) (5) 

has the same norm as the uncorrelated state | $ ) so that the norm in the 
denominator of an expectation value 



does not require special treatment as for example in the Jastrow method [4,8]. 

If the correlator C is chosen to be state independent one may equivalently use 
correlated states | ^ ) and uncorrelated operators B or uncorrelated states 
I $ ) and correlated operators B 

(^1 5 1^') = ($1 ^ 1$') , (7) 
where the correlated operator 

B ^C"^ B C ^C-^ B C (8) 

does not dependent on the particular states | $ ) and | ). 

If C is not unitary as for example in the Jastrow method [4] one may from eq. 
(6) also define a correlated operator Sjastrow ^ 

B _ ^Jastrow C'jastrow /qn 

^Jastrow - I \ r<\ r I tb \ ' ^ ' 

\ ^ I Jastrow^ Jastrow | ^ ) 

but this depends on the state | $ ) . 

For a unitary correlator C there is the duality that one may regard C as a part 
of the trial state C | $ ) which provides further variational degrees of freedom 
beyond those already contained in | $ ) . Or one may view C as a unitary 
transformation of the Hamiltonian which amounts to a pre-diagonalization 

\ H \ = \ H C \ = {^k\H \ . (10) 
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This means that the transformed Hamihonian H = C ^ H C is much better 
represented in the product basis | ) than the original Hamiltonian. 

Both ways of viewing the correlator C are mathematically equivalent. The 
first, where C is part of the trial state, is more suited for the Ritz variational 
principle where one wants to find only the ground state and its energy. The 
second is better adopted to the language of multi-configuration mixing and 
effective operators. In any case the physics contained in C is the same, namely 
to describe the short range correlations. 



2.1 Decomposition of correlated operators into n-body operators 



In order to ensure unitarity the correlator C is written as 



C = expl-iG} , G = G'^ (11) 

where G is the hermitean generator of the correlations. Because we treat 
fermion systems in this work, particle number conserving generators are con- 
sidered only. Furthermore the generator G is invariant under permutations 

of the particles. G has to be a two-body operator or higher because a one- 
body operator would only cause a unitary transformation of the single-particle 
states and this variational degree of freedom is already present in the product- 
state 1$). 



A 

G^Yl j) + three-body + higher (12) 

i<j 

Any correlated operator 

B^C'^ BC^ exp{+iG} B ex.p{-iG} (13) 

will therefore be in general a superposition of zero-, one-, two-, three- and 
higher-particle operators. 

In order to avoid confusion we will use the following notation. An irreducible 
n-body operator B in an 74-particle space is denoted as 



BT = 



A 

E 

ii<---<in 



(14) 
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where • • • , in) is the actual operator with the n particle indices ii,--- , 

Ifn > A then B^^^ = 0. For example, the Hamiltonian which consists of a one- 
body kinetic energy T and a two-body potential V is written in this elaborate 

way as 



A A 



H^T + V^HA^T^2^+vi'^^Y: *™ (0 + E j) ■ (15) 

i i<j 

A missing superscript [n] means that the operator is a combination of different 
many-body operators. Of course later on when there is no danger of confusion, 
superscripts and subscripts are omitted again. 

The decomposition of the correlated operator B = Ba into a sum of irreducible 
n-body operators is written as 



A 

Ba = Ca'BaCa = J2B^a , (16) 

n=0 

where in the A-body space the n-body operator 



E ^'"k^i,---,in) for n>l (17) 

ii<---<i„ 

is given by 



n— 1 

6M:=5n-^5[fl (18) 

k=l 

n— 1 n 

= C-'B^Cr,-Y: E b^'\^l,■■■ ,^k) . 

k=l ii<---<ik 

The first non-vanishing fel"' depends on how many particles the operator B 
connects. The trivial case is when B is proportional to the unit operator then 
Ba = Ca'BaCa = = Ba. H B = E^<...<.^ ^l'"' (^1, ■ ■ ■ , U is an m-body 
operator then the first non- vanishing contribution comes from fc'"*! and is given 

by 

feH = Cm . (19) 

The decomposition (16) is often called "cluster expansion". This is at this 
stage somewhat misleading as the particles ii,i2,-'' whose coordinates 
show up in ^["■^(ii, • • • , i„) need not be close to each other. 
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From now on wc assume the generator to be a two-body operator G — 
J2t<j g{h j)- Of particular interest is the effect of the correlation on one- and 
two-body operators. For example the kinetic energy, a one-body operator, 
transforms to 



fA = C],TACA = f]l^+f^I^+f^I^ + --- , (20) 
where the one-body part is just the regular uncorrelated kinetic energy 



n^-E^H^) = Et^%) = n^ (21) 

i i 

because the generator is a two-body operator. The two-body part T4' is 
obtained from the general eqs. (17) and (18) as 



Tfl = ^?2l(i,j) with (22) 

i<j 

= c{i,j)-'{tm(i)+t^n^j)) c(i,j) - (tW(i) +tW(j)) 

and c(i,j) = exp{— ^(^[^'(i, j)} denotes the correlator between particle i and 
j. The obvious meaning is that for each pair the two-body part is the 

difference between correlated and uncorrelated kinetic energy for the pair. The 
three-body part of the correlated kinetic energy operator is according to eqs. 
(17) and (18) given by 

f|^= Yl t^'\hJ,k) and (23) 

i<j<k 

P\z,j,k) = c{z,j,k)-' (tW(^) +tW(j) +tW(^)) c{i,j,k) 

-(tW(i)+tW(i)+tW(fc)) 

with c{i,j, k) = exp{— i(5f[^l(i, j) + g^'^^(i, k) + g^'^^j, k))}. The three-body part 
correlates all triples {i,j,k) of the system. But it is important to note that 
only the genuine irreducible three-body correlations appear in (23) because 
all two-body correlations and the uncorrelated kinetic energy are subtracted 
out. 

When transforming a two-body operator V^'^ = J2i<j v^'^Kh j) ^i^^ poten- 
tial, the correlated operator 
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Va = Va Ca = t^f +Vf + --- (24) 
starts with a two-body part which is aheady correlated 

v'i^-jlv^'Khj) with (25) 

v^'^\i,j) = c(i, j)"^ v^'^\i,j) c(i,j) . 

This means that the bare interaction is in lowest order replaced by the 
correlated interaction v'^l which will be much less repulsive at short distances. 

The specific form of eq. (18) for the three-body part is 



E v^'\ij,k) with (26) 

i<j<k 

-(y^^\l^j) +v^^\^^k) +v^'\j,k)] 



From the exphcit way of writing eqs. (23) and (26) it is obvious that the three- 
body part will only contribute appreciably if the probability of finding three 
particles simultaneously in the correlation volume is high. The importance of 
the three-body correlations will be discussed in more detail in section 3.3. 

Another very important issue is the "cluster decomposition property". If a 
system decomposes into locally separated subsystems with particle numbers 
A,B,...,Z, which means that the uncorrected state can be written as 

|$)= |$z) (27) 

where | $a ), \^b) and so on are locally disconnected, then 

C\^)^Ca\^a)®Cb\^b)®---®Cz\^z)- (28) 
Prom the construction of C, which is 



C = exp -z 9^"%3)\, (29) 



the validity of the decomposition property is obvious because first, g^^\iij) = 
if i and j do not belong to the same subsystem and second, [g^'^^{i, j), g^'^^k, /)] — 
if (i, j) e A and {k, I) ^ A. Therefore, C can be decomposed as 
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C|$) = exp 1$^)®... (30) 

■■■®expl -iY^g'^^\i,j)\ |$z) 

and hence the "cluster decomposition property" is fulfilled. 

It is easy to see that the cluster decomposition property of the correlator C is 
carried over to each irreducible n-body term in B (see eq. (16)). 

2.2 Representation as coordinate transformation 

The repulsion at short distances will keep the particles from getting too close to 
each other. This feature is not included in a (symmetrized/antisymmetrized) 
product state where the probability to find two particles (with different spin) 
at the same place is just the product of the one-body densities (± exchange 
term). Therefore, the task of the correlator is to shift two particles away from 
each other whenever they are close. 

In this section we define a generator g{l,2) which is producing an outward 
radial shift depending on the uncorrelated distance of the particles. Wc express 
the action of the correlator in the two-body space c(l, 2) = exp{—ig{l, 2)} in 
terms of a coordinate transformation R±{r) of the radial distance r. 

The relative distance and the relative momentum for equal mass particles are 
denoted by 



f = 5(1) - x{2) , g = ^ (p(l) - p(2)) , (31) 

respectively. A generator 5'(1,2) = g{r,q) which creates a position-dependent 
shift may be written in the hermitean form 




r=\f\ . (32) 

One expects the unitary operator c(l,2) = exp{—ig(f,p)} to shift a relative 
position r to about r + s{r)^. The exact transformation will be derived below. 
The idea is to find a suited function s{r) such that g{r, p) moves the probability 
amplitude out of the classically forbidden region of the repulsion. 

In coordinate representation the action of the generator g{r,q) on a relative 
wave function ( r | ) is given by 
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(1 (9 s s d \ 

1 , . d 



^-i^= s{r)— r^s(r) (r|(/)) . (34) 
rJs{r) ' 

Prom that we obtain the coordinate representation of the correlator 



(r I c I 0)=exp <j ^^=s{r)-^ r^/s(r)^ (^10) (35) 



1 . . d 



exp<{-s(r)— > r^s(r) (r|0) . (36) 



The last line can be verified using the power expansion of the exponential. 
Writing c in this way suggests a coordinate transformation 



or ay 

such that c amounts to a shift of —1 in the coordinate y. The transformations 
between y and r are denoted by Y{r) and R{y) 



r^y = Y{r) (38) 
y^r = R{y). (39) 

They are by construction the inverse of each other 



R{Y{r)) = r and Y{R{y)) = y . (40) 
From (37) follow the differential equations 



^y(r) = ^ and R{y) ^ s{R{y)) (41) 

and the integral equation 



rR{y) d£ 

W)' '''' 

which defines the transformation R{y). Rewriting eq. (36) in terms of y and 
r = R{y) one obtains 
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{R{y)-\c\ct>\ 

r 



R{y)^s{R{y)) \ dy\ ' r 

(43) 



The action of exp{ — ^ } is to shift y by —1 and to leave the direction 
unchanged. Rewriting (43) in terms of r yields 



(,|e|,).^(IM^.£WlM(fi(yM-l)I|« (44, 
r \ s{r) r 

or for the inverse transformation 



(r I c \ — 



R(Yir)^,) MMYirHD) ^niY(r)^,^f 

r \^ s{r) r 



These two equations show that the transformations R and Y appear only in 
the form R(Y{r) ± 1). Therefore it is convenient to introduce the functions 



i?±(r) = R{Y{r)±l) 
with which we also express the terms under the square root as 



(46) 



or s[r) 
Prom the integral equation (42) follows 



(47) 



«(2/±l) d{ 



R{y) s(0 



fR± d^ 

y±l-y or / — — = ±1 
Jr s(0 



(48) 



which is the formal definition of R±{r). From the fact that Y{r) and R{y) are 
inverse of each other follows 



R±{R^{r))=r , 
which means that R+ is the inverse of R- and hence 



(49) 



it:4(r) = \R'{R±(r)) 



-1 



(50) 
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With help of the coordinate transformations -R+(r) and R-{r) the correlated 
wave function for the relative motion of two particles given in (44) can be 
written as 



{r\c\(f)] 



R_{r) 



(51) 



and the inverse transformation as 



(r|c 



(52) 



Hence, the unitary correlator c = exp{— igf} is uniquely given by the coordinate 
transformation i?_|_(r). 

With the correlated state given in (51) one can define the correlation volume 
Vc by the square of the defect wave-function which is the difference between 
a uniform uncorrelated state ( r | ) and (r|'0) = (r|c|0) 



K := 47r P dr r- R^{r)^J RL{r) 

i 



(53) 



47r / dr 



t-R^{t)JR'+{t) 



A typical case is illustrated in fig. 2. 

By looking at a matrix element {ipi \ b{r) \ip2) of a local operator b(f), which 
acts on the relative distance, the meaning of the different parts in (51) and 
(52) becomes clear 



( -01 1 b{f) I ■02 ) = j r^dr dil ( -01 1 r ) b{f) ( r | -02 
= J r^dr dQ {<pi\ | ^ ) b{r) ( r | c | 02 ) 



(54) 



rMrdn 



R(r) 



R'_(r) (V'i|i?-(r)-)6(r) ( i?_(r)- | ^^2 ! 



J r!dr_dQ_ {^i\f-)b f i?+(r_)^ j ( r- | V2 



The factor ^ ^ RL{r) in (51) is just the square root of the Jacobian of the 
transformation 
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dV_ = rldr_dn_ = (^^^^ R'_{r) r^drdil , (56) 
where dJl = dfl- = sin OdOdcp remains unchanged. 

Thus, the correlated wave function (51) at a relative distance r is given by the 
uncorrelated one taken at a closer distance i?_(r)^ times the square root of 
the Jacobian which takes care of unitarity. 

Equivalently one can use the correlated relative distance 



9 ^ r c ^ R+{r)^ . (57) 

Due to the unitarity property c^c = 1 any local operator b{r) can be readily 
transformed to the correlated one 



b{r) ^ cH{r) c ^ c-^ b{r) c ^ b \^R+{r)^j . (58) 

The correlated operator b{f) is just the uncorrelated one taken at the trans- 
formed relative distance R+{r)^. This is also the physical meaning of the last 
line in eq. (54). 

The correlated radial momentum can be calculated in analogy to (54) by 
using the coordinate transformation (55) r — > r_ and the relation (50) between 
R'_ and R'^. The result is 



Qr = C ^ Qr C = . - Qr f i (59) 

where 



r d 

{f\qr\(f))^{r\--q\(f))^-i-^ {r\(f)) ■ (60) 

The correlated relative angular momentum is the same as the uncorrelated 
one 



T=c ^lc = c ^rxqc = l. (61) 
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r [fm] r [fm] 



Fig. 1. L.h.s.: typical coordinate transformation R+{r) and its inverse 
R-{r); r.h.s: shifts r — R-{r) (full line) and R+ir) — r (dashed line). 

The present correlator conserves angular momentum and spin for fermions. 
The realistic nucleon-nucleon interaction, however, induces strong tensor cor- 
relations. For such a situation, which will be subject of a forthcoming paper, 
the coordinate transformation R-{r, ai, (72) will not only depend on the rela- 
tive distance |r| but also on the angles between r and the spins ai and a2 of 
the two particles. 

An intuitive relation between the coordinate transformation R±{r) and the 
original s{r) in eq. (32) can be obtained from the first term of the Taylor 
expansion around y using (41) 

R^{r) = R{Y{r) ± 1) = R{Y{r)) ± R'{Y{r)) + ■■■ (62) 
= r ± s(r) + ■ • • . 

For s(r) small compared to r, s(r) is the amount by which the two particles 
are shifted away from their original distance r, as anticipated earlier. 

Instead of discussing the shift function s(r) as the variational degree of freedom 
we will consider the equivalent coordinate transformation R±{r) directly. 

In order to avoid unphysical transformations, R±{r) has to increase monoton- 
ically with r and has to be differentiable everywhere. For correlations of finite 
range -R±(r) r for large r. A typical example which will be used later is 
shown in fig. 1 in terms of R±{r) on the l.h.s.. The r.h.s. displays the shifts 
r — R-{r) and R+{r) — r which is the essential information. The shifts rise 
rapidly at r = to a maximum of 0.15 fm around r = 0.4 fm. This results in 
a strong depletion of the relative wave function as can be seen in fig. 2 where 
a uniform uncorrelated state is compared with the correlated one as given in 
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[MeV] 
100 

50 



-50 



0.5 1 1.5 2 2.5 

r [fm] 

Fig. 2. Uncorrelated state {r\(j)) (dashed line), correlated state {r\tp) 
(full line) and potential v{r) (thin line) as a function of the relative 
distance r. 

eq. (51). The corresponding potential is also included in the figure. 

2. 3 One- and two-body part of the correlated Hamiltonian 

Of particular interest is the n-body expansion of the correlated Hamiltonian. 
Starting from eq. (58) the two-body part of the correlated potential is given 
by 

t;Pl(z, j) = c(z, c{t,j) = V (/?+(r,,)) (63) 

where rij = x{i) — x{j) and rij = \rij\. 

The three-body part v^^^{i,j, k) assumes the more complicated form given in 
(26). It will be discussed in section 3.3. 

The correlated kinetic energy contains the uncorrelated kinetic energy as the 
first term 

?^l(z)=t(z) (64) 

because g^'^^{i,j) is a two-body operator and hence Ci = 1 in the one-body 
space. The two-body part can be calculated by first separating the kinetic 
energy of two particles into relative and centre of mass energy 
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where q^j — {p{i)—p{j))/2 is the relative momentum. Since g^'^^{i,j) commutes 
with the centre of mass coordinate, only the relative momentum is correlated 



c(i,j)-'(t(i)+t(j)) c(i,j) (66) 
= - c(i,j)-^ ql c{i,j) + — {p(i) +p{j)f . 

Decomposing qf, into a radial and angular momentum part 



1 -'2 

ql = 4ij<lrij + jihj > (67) 
where qr = ^ -q and I = r xq, one can utilize eqs. (59) and (61) to obtain 

c{i,j)-'qlc{i,j) (68) 



Commuting the two radial momenta qrij in the first term on the r.h.s. to the 



left and to the right of rij/ J R!^_{rij), respectively, yields 



c{i,j)-'q^.c{z,j) (69) 



1 1-2 



(i?;(r,,))2 + (i?+(n,))2"^'^- + '"^'''^''-^'^ ' 



where 



u 



1/m 



2 Kir^, 



5 fRlin 



+ l^^rP4 \ . (70) 



R'^{r) and R"'{r) are the second and the third derivative of R+{r) with respect 
to r. Thus, the two-body part of the correlated kinetic energy may be split 
into three parts 



(71) 



irij 



m 



- 1 



qrij ~l~ 



m 



{R+{rij)Y r\ 



^3 
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The first is a two-body radial kinetic energy with a reduced radial mass 



m 

^^r{r) = — 



-1 



(72) 



which depends on the distance r between the particles in the pair. The second 
part is a relative rotational energy with a reduced angular mass 



m 



n -1 



1 



(73) 



The third part SPJ(i, j) is a local two-body potential which depends on deriva- 
tives of -R+(r) as given by eq. (70). All three parts vanish outside the correla- 
tion range when R+{r) — > r. 



The one- and two-body parts of the correlated Hamiltonian add up to 



^[1] + H^'^] = J2 t{i) 



(74) 



rij 



Qrij ~l~ 



A 

+E 

i<j L 



f 

2 'ij 



The first part is the uncorrelated one-body kinetic energy. The second is a 
momentum dependent two-body interaction which is restricted to the corre- 
lation area. The third is the tamed two-body potential which has no repulsive 
core anymore. Its shape will be determined in the following section where the 
choice of R+{r) is discussed. The taming of v{r) has to be paid by the mo- 
mentum dependent two-body interaction. In Cartesian coordinates this part 
may also be written as 



ql 



where 



1 



qr + 



Mr) 



q + ql 



Mrir) 



(75) 



m 
2" 



1 



(i?;(r))2 {R^{r)y 



(76) 
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2-4 Choice of the correlation function R±{r) 



The correlator is a functional of s{r) or equivalently of the coordinate trans- 
formation i?-(r) 



I*) = C[i?_] 1$) . (77) 

Therefore, R-{r) represents additional variational freedom and according to 
the Ritz variational principle R-{r) could be determined by minimizing the 
ground state energy 



( ^ I iJ I ^ ) = ( $ I C^[R^] H C[R_] I $ ) (78) 

not only with respect to | $ ) but also with respect to R^{r). 

However, since only the first terms in the n-body expansion of the correlated 

Hamiltonian H — W^^+W"^^ H (see eq. (16)) will be used, the exact ground 

state energy is in general no lower bound on ( $ | H^^^ | $ ) + ( $ | | $ ) 
anymore because contributions ( $ | iJ'"^ | $ ) for n > 3 can be positive and 
negative so that only the complete sum up to the total particle number A has 
the lower bound. 



In order to make the three- and higher-body correlation part as small as 
possible one should use the freedom in the separation of short and long range 
correlations, i.e. what should be taken care of by C[i?_] and | $ ) respectively, 
to find an optimal correlated two-body Hamiltonian. One expects the three- 
body part of the correlated Hamiltonian which is induced by C[i?_] to be small 
if the correlation range in which R_{r) ^ r (or s(r) ^ 0) is short compared 
to the mean distance between the particles. A measure of the smallness is the 
uncorrelated density po times the correlation volume Vc, defined in (53). In 
section 3.3 where we calculate three-body terms a relation between poK ^-nd 
the three-body correlations is given. Prom that one can judge how small poK 
actually should be. 

In the following we investigate for a two-body system the general idea that 
the correlator C can be imagined as a unitary mapping of uncorrelated states 
onto correlated states which approximately diagonalize the Hamiltonian H. 



H\7Pe)^E\^e) (79) 

\iPe) are the exact bound and scattering eigcnstates of the Hamiltonian 
H — j^q'^+v{r) for the relative motion. The repulsive core of v{r) will suppress 
the wave function ( r | ■0b ) for small r. The goal is to find for a range of energies 
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an optimal correlator c[R-] which removes these repulsive correlations from 

If one demands that in a range r < A the correlated trial state and the exact 
state coincide 



(r| c[i?_] 10) = (rlVe) for < |r1 < A (80) 
one obtains for angular momentum equal zero from eq. (51) 



./RLir) u{R4r)) ^ UE{r) (81) 
or the integral equation 



where i?_ (r = 0) = is assumed and the radial wave functions are denoted as 

UE{r) — r{f\ipE) and u{r) — r{r\(l)) . (83) 

As already discussed in a more general frame in section 2 the correlator C or 
the coordinate transformation R_{r) as defined in eq. (82) depend on the trial 
function u{r) and on the Hamiltonian via its eigenfunctions UE{r). 

For systems at small energies a first attempt is to assume the trial state 
(r I 0) = 1/\/Nq to be uniform across the repulsive core [11]. For this special 
case the coordinate transformation R_ (r) is given by the differential equation 



«iM^iV.(^)\ (84) 
where the norm Nq is still to be chosen. 

The dependence of i?_(r) on A^o; which is displayed in fig. 3 l.h.s. in terms of 
the deviation r — R-{r), will now be used to limit the range of the correlation. 
For this the MTV potential ([13] and eq. (A.l)), already shown in fig. 1, is 
chosen as an example and the energy E is set to zero. Among the set of curves 
there is one which touches the a;~axis at a distance A which means 



i?_(A) = A and R'_{X) = 1 . (85) 
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Fig. 3. L.h.s.: Coordinate transformation R-{r) plotted as r — R-{r) 
for E = and different Nq. R.h.s: Exact radial wave function UE=oir) 
(solid line) and trial function u{r) for which i?_(A) = A and R!_{^) = 1- 

From the Schrodinger equation for angular momentum zero 
m (f (r) — E) Ue{t) 



one sees that at the classical turning point v{rt) = E the curvature of UE{r) 
switches sign so that a linear radial function u{r) = t/^/Nq can be found 
which intersects the exact solution twice, see fig. 3 r.h.s., such that at r = A 
trial and exact state coincide 



u{\) = ue{X) 
and their norms up to A are also equal 



^7) 



A nX 

dr {u{r)f = / dr {uE{r)f . (88) 
Jo 

The transformation R-{r) which fulfills the conditions (85) or equivalently 
(87) and (88) provides a natural division into a short and a long range part. 
This division is an inherent property of all methods which derive effective 
interactions [12]. At r = A one may split c[i?_] into 



c[R_] = c[R{_] c[R'J] 
where 



^9) 
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, i?_(r) for < r < A „ r for < r < A 

Ri = { ^ ^ - - and RH = 1 ^ ~ (90) 

r for A < r I R-{r) for A < r 

because c[R{_] commutes with c[J?:L^] if (85) is fulfilled. For r > A the uncorre- 
lated state ( r | ) is supposed to describe the long range part so that we set 

The schematic case of a uniform uncorrelated state illustrates nicely how the 
method works. As the correlation function i?_(r) is chosen such that the cor- 
related state equals the exact solution for a given energy E (see eq. (80)) it 
also solves the Schrodinger equation (86) 

(^1 (^^Ql(lr + v{r)^c[R-] \(p)^E{f\ c[R-] |0) (91) 



or 



(r| (—c[R-]-^qlqrc[R-] + c[R-]-'^ v{r) c[R-]) \ (p) = E {r\ (p) .{92) 

Since (r| 0) = I/^/Nq is uniform, the relative kinetic energy reduces to (c.f. 
eqs. (69) and (70)) 



{f\^c[R.]-'qlqrc[R.] |0)=SP](r)-^ (93) 
and hence 



SPl(r) + ^^^^V) for r < A . (94) 

This means that a correlator obtained from eq. (84) for £■ = transforms 
the original Hamiltonian for small distances r < A, where the potential is 
strongly repulsive, into a correlated Hamiltonian with a vanishing local po- 
tential M[^l(r) +{;[^l(r) = 0. Thus the repulsive core is completely absent in the 
correlated Hamiltonian. 

If the potential is purely repulsive {v{r) > for all r) there is no natural 
scale for the division into long and short ranged components. In this case the 
quest for a correlation volume small enough for the two-body approximation 
will determine the range of the correlation. Further discussions on repulsive 
potentials can be found in section 4.1 where we treat spin-isospin-dependent 
correlators. 

The next point is the dependence on the energy. Fig. 4 displays r — i?_(r) 
for a separable correlation at different energies between and SOMeV which 
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Fig. 4. Energy dependence of the coordinate transformation R-(r), 
E = (full line), E = 40MeV (dashed line), E = 80MeV (dotted 
line), rt classical turning points, A range of correlators. 



is the kinetic energy of two nucleons colliding with opposite Fermi momenta. 
As the uncorrelated state we choose the spherical Bessel function so that 
the radial function u{r) = rjo( V mEr) / \/Nq = sin( V mEr) / ^/Nq is the I = 
solution of the free Schrodinger equation. It turns out that for distances r < rt 
the correlation function i?_(r) is rather energy independent. For r > rt one 
expects an influence from the uncorrelated state so that the shape of R-{r) 
may differ for vt < r < X from the one shown in fig. 4 depending on the choice 
of the uncorrelated state. 

The kink in the wave function at r = A (see fig. 3 r.h.s.) caused by cutting 
off i?_(r) at r = A in a non-analytic way should be removed. Therefore, the 
short range part R^_{r) will be parameterized by an analytic function which 
approximates R-{r) for < r < vt rather well and leaves some freedom for 
r > rt for adjusting to the actual trial states used. 

Since later on in the many-body case the product state will consist of Gaus- 
sian shaped single-particle states, the uncorrelated state of the two-particle 
system, from which the coordinate transformation i?_(r) is deduced, is also 
assumed to be a Gaussian at small r (and not a constant or a Bessel function) 
which joins smoothly to an exponential radial function at r = p 



exp|— ^^r^l for r<p 
(r|</>) = ^_ '^'^ \ -V (95) 

'Nip) [£exp{^-^} for r>p 

The length k = (— mi?o)^^^^ is chosen such that the long range part, where 
the correlator should not modify the trial state, has the correct fall off. 
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Fig. 5. L.h.s.: MTV potential v{r) in MeV, uncorrelated wave func- 
tions for p = 2.62 fm (thin dashed line) and p = 2.5 fm (thick dashed 
line) and corresponding correlated ones which both equal the exact one 
(thick solid line). R.h.s: correlation functions in terms of r — i?_(r) for 
uncorrelated states which map on the exact ground state. 



Fig. 5 displays on the left hand side the uncorrelated trial state ( r | ) for 
p = 2.62 fm and the exact ground state (r| together with the MTV 
potential. The exact solution is pushed out of the repulsive region compared 
to the uncorrelated state. The coordinate transformation i?_(r) obtained from 
(82), which transforms {r\(f)) into (r | ipE ), is shown on the r.h.s. in terms of 
the shift distance r — R_{r) (thin line). Due to the attraction this shift turns 
out to be negative beyond r = 1.2 fm so that more probability is accumulated 
around r ^ 1 fm where the potential is most attractive. This long range effect 
should not be corrected by the correlator but an improved trial state should 
be employed, otherwise the three- and higher-body terms in the correlated 
Hamiltonian will be too large. 

Similar to the uniform trial state discussed above we vary the norm N{p) by 
moving the matching point p until we achieve a separable correlator which is 
the case for p = 2.5 fm (thick line in fig. 5 r.h.s.). This correlation function 
differs only little from the one where trial and exact state coincide for very 
large r {p = 2.62 fm). The resulting shifts r — R-{r) are displayed on the r.h.s. 
of fig. 5. For r < Tt = 0.6 fm they are identical and similar to the case of 
a uniform uncorrelated state. The separation point A turns out to be inside 
r = p, so that the correlation acts only on the Gaussian part. 
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3 Few— body systems with spin— isospin— independent forces 



In this section we apply the concept of the unitary correlation operator to the 
two-, three- and four-body systems and ^'^He. The correlation function 
R+{r) is first parameterized and then the binding energies of the three nuclei 
are minimized with respect to the parameters of i?+(r). It turns out that the 
resulting correlation functions are almost identical. This shows that the general 
concept of a state-independent correlator discussed in the previous section is 
applicable. The last subsection investigates an approximate treatment of the 
three-body term in the correlated Hamiltonian. 

Throughout this section the Malfiiet-Tjon V potential (MTV) [13] is used, 
which is a central interaction without spin and isospin dependence (see ap- 
pendix eq. (A.l)). This potential is widely used in the literature, for instance 
to compare different theoretical many-body approaches [14,15] because it has 
a repulsive core and an attractive long range part chosen by considering the 
phase shifts of the scattering problem. However, already for ^He it overbinds. 
For larger nuclei it gives completely unreahstic results. Therefore another po- 
tential will be used in section 4 where nuclei up to A = 48 are investigated. 



3.1 Parameterization of R±{r) 



For the two-body system everything can be calculated with a correlation func- 
tion R+{r) that is derived as described before and at least numerically known 
at all r. But for many-body systems it is advisable to have an analytic form 
with a limited number of parameters. Therefore we parameterize the numeri- 
cally determined R+{r) by 





R+{r) = ^ + « exp< — exp H • (96) 



The parameter a controls the overall amount of the shift and /3 the length 
scale. T) determines the steepness around r = 0. The double-exponential fall 
off is required to reproduce the exact correlation function in the vicinity of 
the separation distance r = A, which has been derived in the previous section. 
Figure 6 compares the exact correlation function for the bound state with 
the parameterization (96), where a — 0.935 fm, P — 0.95 fm and 77 = 0.39. A 
similar fit can be achieved with a Gaussian replacing the double-exponential 
fall off. 



26 



[fm] 



0.1 



0.05 





0.5 1 1.5 2 2.5 

r [fm] 

Fig. 6. Exact correlation function R+{r) for the ground state (thin line) 
and the parameterization (96) plotted as R+{r) — r. 

3.2 '^H and ^^^He with the Malfliet-Tjon V potential 

In this section the particle number dependence of the correlation function is 
investigated. For that the three parameters a, (3,1] of the ansatz for R+{r) 
given in eq. (96) are varied in order to minimize the binding energy of the 
two-, three- and four-nucleon system. The uncorrelated trial states are kept 
as simple as possible because the short range correlations are supposed to 
be described by the unitary correlator alone. Therefore, the uncorrelated ^He 
and ^He states are taken to be a product of three and four identical Gaussians 
in coordinate space, respectively. The antisymmetry is taken care of by an 
appropriate spin-isospin part of the many-body state. This excludes for ^He 
the admixture of relative p-states which are important for a correct ground 
state energy. 

Thus all three uncorrelated many-body wave functions are characterized by 
only one variational parameter, which for is p as defined in (95) and for 
^'■^He is the width a of the Gaussian shaped single-particle state. 

The correlated Hamiltonian is calculated up to its two-body part as defined 
in section 2.1 

H^^ = Tint + r'^' + V^^^ (97) 
where the intrinsic kinetic energy 
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T,nt = T-T,m, '^-rn = ^I^J2my (98) 

is the kinetic energy minus the centre of mass kinetic energy Tcm- 

Table 1 summarizes the result of minimizing the expectation value of the 
correlated intrinsic Hamiltonian H^^. The first and very important result is 
that the correlation functions R+{r) (represented by their parameters a, P, rj), 
which result from minimizing the energy, are very similar for the two-, three- 
and four-body system. This can also be seen in fig. 7 where they are plotted 
together with the -R+(r) which results from mapping the uncorrelated deutcron 
ground state (95) to the exact relative wave function up to r = A as defined 
in (82) and discussed in section 2.4. 

Since all correlation functions turn out to be almost identical, the initial con- 
cept of a unitary state-independent correlator, which does not depend on the 
actual system it is applied to, is strongly supported. 

The second result is that despite of the simplistic nature of the uncorrelated 
states the binding energies of the A = 2 and A = A systems are very close to the 
exact values. The A = 3 system is still 1.8 MeV above the value of a Green's 
function Monte Carlo calculation [16]. Responsible for that are the missing 
exponential tails in the Gaussian radial functions and a missing component 
in the uncorrelated state in which the spins of the two protons are coupled to 
one and the relative motion is in a p-state. 

One should keep in mind that the total energy results from the subtraction 
of two large numbers, the kinetic and the potential energy, both depending 
sensitively on i?_|_(r). As will be demonstrated with an example below, the 
potential energy almost doubles when including the correlator. 



A 


a 




V 


a 


P 








exact 




[fm] 


[fm] 




[fm2] 


[fm] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


2 


0.939 


0.996 


0.394 




2.42 


4.78 


-5.17 


-0.40 


-0.41 


3 


0.944 


1.018 


0.389 


2.06 




34.97 


-41.49 


-6.52 


-8.26±0.01 


4 


0.954 


1.121 


0.369 


1.51 




76.37 


-107.2 


-30.78 


-31.3±0.2 



Table 1 

Parameters o:,l3,r] of the correlator and a,p of the uncorrelated many- 
body state which minimize the energy ( H'^'^ ) . { Tj^^ ) = { Tint + T^'^^ ) 
is the intrinsic correlated kinetic energy and ( V^"^ ) = {V^'^^) the 
correlated potential energy. The exact binding energies are taken from 
[16]. 
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Fig. 7. Comparison between exact correlation function for the ground 
state of (thick line) and parameterized functions obtained by min- 
imizing the ground state energies of (thin line), ^He (long dashed 
line) and ^He (short dashed line) in terms of R+{r) — r. 

Since the difference between the correlation functions is so small a standard 
correlator with the parameters 

a = 0.94fm /5 = lfm, t] = 0.37 (99) 

is used for all three systems and the minimization is done only with respect 
to a single parameter, p or a, which determines the radial wave functions. 
As expected their values differ very little from those given in table 1. The 
resulting energies are decomposed into their components in fig. 8. 

For each system the uncorrelated expectation values of kinetic, potential and 
total energy are compared with their respective correlated counter parts. The 
uncorrelated energies are calculated with the uncorrelated states which min- 
imize the correlated energy. In all three cases the kinetic correlation energy 
( TI^I ) (indicated by a vertical line) is much less than the potential correlation 
energy ( ) — {V) so that the sum of both reduces the uncorrelated total 
energy from positive values (unbound) down to almost the exact ground state 
energies. 

Figure 9 illustrates the action of the unitary correlator. On the left hand 
side one sees that the repulsive maximum of r^f (r) is reduced by almost a 
factor of three, while the attractive region is enlarged. This corresponds to the 
shift of the wave function away from the repulsive core towards the attractive 
region, c.f. (54). The correlated kinetic energy contributes a potential like 
part Mt^l(r) which is given in eq. (70) and depicted in the middle of fig. 9. For 
short distances r<l fm it has a pattern opposite to the correlated potential. 
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Fig. 8. Energies for A = 2, 3, 4 and MTV potential. Left columns: expec- 
tation values of kinetic energy ( Tint ) , potential ( V ) and sum of both 
{H) . Right columns: expectation values of correlated kinetic energy 
( Tint + r'^^ ) , correlated potential ( yt^] ) and sum of both ( H*^^ ) . 
Arrows indicate exact values from [16]. 
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Fig. 9. L.h.s.: correlated (^[^^(r) = v{R^{r)), full line) and uncorre- 
lated {v{r), thin line) MTV potential. Centre: potential like part of the 
correlated kinetic energy u^'^\r). R.h.s.: MTV potential (thin line) and 
"tamed" potential u[^](r) + wPl(r). All potentials are multiplied with r^. 



attractive at short and repulsive at larger values of r, so that the sum of both 
^[2](^) _|_ {;[2](^) (r.h.s. of fig. 9) almost vanishes in the region of the repulsive 
core of the MTV potential. As discussed in section 2.4 a uniform uncorrelated 
state which is mapped by the unitary correlator c onto the exact eigenstate 
of the Hamiltonian at energy E leads to u^'^^{r) + v^'^^{r) = E for small r. This 
explains the small values of vP^{r) +v^'^\r) at small distances. The deviations 
from a constant E are due to the fact that R+{r) is parameterized (96) and 
the trial state is not uniform. 
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Fig. 10. L.h.s.: Inverse reduced mass in radial direction. R.h.s.: Inverse 
relative moment of inertia. The values are multiplied with and in 
units of m/2. 

The r.h.s. of fig. 9 shows that the repulsive core of the original potential 
is completely transformed away by the unitary correlator. The uncorrelated 
states feel only a tamed potential which is purely attractive. 

Since the correlator creates "holes" in the many-body state it induces the 
additional kinetic energy One part of that is u^'^\r), the other part can be 
understood as a momentum dependent two-body potential which is written 
in eq. (69) in terms of a reduced mass and a reduced moment of inertia. Their 
dependence on the relative distance r is displayed in fig. 10. The l.h.s. shows 
the inverse of the reduced mass for the radial direction as calculated from 
eq. (72) in units of m/2 and multiplied by r^. It results in an increase of 
the correlated kinetic energy. The r.h.s. of fig. 10 depicts the inverse of the 
reduced moment of inertia times which is given by eq. (73). For relative 
angular momenta / = this part does not contribute. But for components 
with Z 7^ it is always negative. This is due to the fact that the correlator 
shifts the particles away from each other so that their moment of inertia 
increases or the inverse of that decreases. Therefore, the correlation energy 
which is the difference between correlated and uncorrelated energy has to be 
negative. Altogether the expectation value of the momentum-dependent part 



of f = Y.iKj^'^Khj) is for A = 4 only 3.3 MeV compared to 12 MeV for 



3.3 Three-body part of the correlated Hamiltonian 

The ansatz for the correlator discussed in section 2 includes only a two-body 
generator g^'^^(i,j). Nevertheless the correlated operators contain three- and 




31 



higher-body operators because g^^\i,j) appears in the exponent (see eq. (26)). 
In principle it is possible that the nuclear many-body system requires also 
genuine three-body correlations so that the corresponding unitary correlator 
would be written as 



C = exp(-^^5p](i,j)-^ E 9^'\hj,k)) . (100) 

^ i<j i<j<k ' 

In the n-body expansion of a correlated operator B = C~^BC the three- 
body correlations induced by (jf''^] would appear first in the three-body part 
together with the three-body part originating from g^^^ It will be very hard 
or even impossible to distinguish between the two contributions just from 
the three- or four-particle ground state energies. This uncertainty together 
with the complicated calculation of the three-body contribution of (/[^' as 
given in (26) suggests to estimate the three-body contribution first before 
doing elaborate calculations. A small three-body contribution to C^^BC from 
C = exp{— i X]^^- 5ft^l(i, j)} is of course desirable for justifying the two-body 
approximation. On the other hand if it turns out to be small even a small 
genuine three-body generator ^fl^l might be of equal or more significance and 
one would have to discuss its shape and importance. Thus the idea is that 
either it is sufficient to consider the correlated operators only up to the two- 
body part or the method becomes too complicated for practical purposes. 

The two-body approximation will be tested by calculating the correlated 
Hamiltonian up to three-body operators with the following approximation 
which is a partial sum of a power expansion of C~^HC [17]. 

c{i,j,k)-' hit,j,k)c{i,j,k) (101) 
~ 4 ^ ^) ^('^^ 

P{a,l3,l} 

where h{i,j,k) = t^^^(i) +t^^^{j) +t^^^{k) +v^^\i,j) +v^^^{i,k) +v^'^^{j,k). The 
greek indices Q;,/3, 7 denote the three possibilities to pick a pair of particle 
indices out of {i,j,k}. The summation runs over all 3! possibilities to or- 
der the triple {q;,/3, 7}. c{a) is the two-body correlation operator acting on 
the pair labeled by a. This approximation correlates all pairs — so to speak 
— one after the other in all different orderings instead of simultaneously as 
exp{-i Etcj g^'^Kh j)} would do. 

It is not exactly unitary anymore but it conserves the norm in the sense that 
the correlated unit operator remains unchanged 



1^3 = 1 . (102) 
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The actual calculations for the three-body and four-body system are done by 
Monte Carlo integration. The derivatives of the kinetic energy are calculated 
numerically. 

For example in the three-body system the potential energy is given by 

{y''')=h E f d'xid'x2d'x3 (103) 

3 

X v{xi - Xj) I 52,53 I c{a) c{l3) c(7) | $)|^ 
and the kinetic energy is calculated as 

{T''')-l E U'xA'x,d?x, (104) 
X ^ E I ( ^1' ^2, ^3 I c{a) c{(5) c(7) I $ ) r ■ 

TTl j Xi 

The uncorrelated states | $ ) are those which minimize the two-body cor- 
related Hamiltonian together with the standard correlation function R+{r) 
defined in section 3.2. 

The results are summarized and compared with the two-body approximation 
(for which analytical expressions exist) in table 2. 

The correlated three-body kinetic energy ( ) is for both nuclei within the 
error bars of the Monte Carlo integration the same as the correlated two- 
body kinetic energy ( T^^ ) . This means that the three-body contribution 
( f [3] ) = {fgl) - { fgl ) is within ±0.06 MeV for A = A. 

The three-body contribution to the correlated potential ( V^^^ ) — { V*^^ ) — 
{ ) is of the order of 0.14 MeV for A = 3 and 1.09 MeV for A = A. Both 
corrections are small compared to the difference ( V'""^ ) ~ ( ^ ) between the 
two-body correlated and the uncorrelated potential which is —18.08 MeV for 
A = 3 and -55.89 MeV for A = A. 

It is interesting to note that the three-body contribution shifts up the energy 
again slightly by 2% of the gain in binding achieved by the two-body part of 
the correlation. However, the uncorrelated state is not varied for the three- 
body approximation separately but taken from the energy minimization within 
the two-body approximation. A variation could lower the energy again. 

The three-body part of the correlated potential is positive in all cases we con- 
sidered. The reason is that a simultaneous action of c{i,j, k) on three particles 
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which are close shifts them further away from each other than the sum of pair 
wise shifts c{i,j), c{i, k) and c(j, A;). As c{i. j) is optimized to move probabihty 
from the repulsive into the attractive region (see fig. 5) a further shift results 
in a loss of binding. 

In order to see the dependence on density of the three-body contribution to the 
total kinetic energy and potential energy we calculate (T^^^) and (V^^^) as a 
function of the width parameter a for ^He. In fig. 11 the ratios ( T'^' ) / ( T^^ ) 
and ( V^^^ ) / ( V'-'^ ) between the three-body part and the total correlated 
energies are plotted versus k = PmaxVc, where pmax is the maximum of the 
one-body density and Vc = 0.16fm"^ is the correlation volume defined in eq. 
(53) for the correlation function (parameters in eq. (99)) used in this section. 
While the three-body part of the correlated kinetic energy stays rather small 
the three-body part of the correlated potential energy grows linearly with k. 
The ratios stay below 1% for < k < 0.1. The ground state values for ^He 





^ = 3 


Error 


A = A 


Error 


a [fm^] 


2.07 


1.55 


( Tint ) 


30.05 




60.39 




/T'C'2 \ 
\-^int 1 


35.12 




75.58 




K-'-int 1 


35.13 


±0.02 


75.56 


±0.06 


{ f PI ) 


5.07 




15.19 




( Tt^l ) 


0.01 


±0.02 


-0.02 


±0.06 


(V) 


-23.54 




-50.23 






-41.62 




-106.12 






-41.48 


±0.001 


-105.03 


±0.002 




-18.08 




-55.89 






0.14 


±0.001 


1.09 


±0.002 


{H) 


6.51 




10.16 






-6.50 




-30.54 






-6.35 


±0.02 


-29.47 


±0.06 




-13.01 




-40.70 






0.15 


±0.02 


1.07 


±0.06 



Table 2 

Uncorrelated energies, two- and three-body approximation of corre- 
lated energies, all in MeV. 
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-0.004 



-0.008 



0.02 0.04 0.06 0.08 
/c 

Fig. 11. Ratios of correlated three-body part to total kinetic 
( f [3] |> / ^ fC3 ^ potential energy ( f^t^] ^ j ^ yCZ ^ ^ function of 

1^ = PmaxVc {Vc = O.lGfm^) for ^He. Values of k for the ground states 
of ^He and ^He are indicated. Lines are only guiding the eye. 

and ^He are in this range. One should however keep in mind that the error 
induced in the total energy ( T"-^^ ) + ( V'-^'^ ) by neglecting the three-body 
contributions is larger because kinetic and potential energy have opposite sign 
and are of similar size. A rule of thumb is to stay with k below 0.1 if an 
accuracy of less then 1% is desired for the effective potential in two-body 
approximation. 

Since the estimated three-body contribution turns out to be so small com- 
pared to the leading two-body contribution it is not worthwhile to pursue it 
further. The other uncertainties, like a genuine three-body contribution, the 
parameterization of the correlation function R+{r) or the simplicity of the 
uncorrelated many-body state are too large. 



4 Many— body systems with spin— isospin— dependent forces 

In this section the unitary correlation operator method is applied to nuclei up 
to mass A = 48. Since the Malfiiet-Tjon V potential overbinds nuclei beyond 
A = A strongly the so called modified Afnan-Tang S3 potential [18,19] for 
which several calculations exist in the literature [19-22] is used. It is specified 
in the appendix in eq. (A. 2). The radial dependence of the potential in the 
two even-parity states with S = and S = 1 and in the odd state is depicted 
in fig. 12. 

The uncorrelated many-body states are Slater determinants composed of 
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single-particle states 



( X I gfe ) = exp ^ — — \ I ms,k ) | mt,k ) (105) 

with a Gaussian-shaped coordinate space part, a two-component spinor 
I rris^fc = ±1 ) and an isospin part | rrit^k = {proton, neutron} ). The complex 
parameters bk and are the variational degrees of freedom of the uncorrelated 
many-body state | $ ). They contain the mean positions, the mean momenta 
and the complex width of the packets. For details see ref . [23] . 

All proton and neutron spins are paired to zero except if the number is odd. 
4-1 Spin-isospin-dependent correlator 



Analogue to the central potential the generator of the unitary correlator is 
decomposed into the four spin-isospin channels 



{0,1} 

^Pl(r,g,ai,a2,Ti,T2) = ^ 9skr,q) Hs^^t (106) 

S,T 

with the projection operators Ils=o = |(1 — <?i<?2), n5=i = ;|(3 + 171^2) and 
11^=0, nT=i accordingly. 

The unitary correlator in A-hody space is herewith 



C^ = expi - iJ29^'^\'i'khQki^^k,^i,rk,ri)> . (107) 

k<l 

For the two-body approximation of the correlated operator only c = C2 in the 
two-body space is needed 



exp - i ^ ^§(r, q) Us^uA (108) 

S,T ^ 



{0,1} 

= Yl e^p{ - WsW, g) } Hr . 

S,T 

Since Us and are projection operators the correlator becomes a direct sum 
of four commuting unitary correlators, one for each channel. Compared to the 
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U , I , I , I , I , I , U U , 1 , 1 , 1 , 1 , 1 , U 

1 2 3 1 2 3 

r [fm] r [fm] 



Fig. 12. L.h.s.: Radial dependence of the potentials for even parity chan- 
nels (5, T) = (1, 0) and (0, 1), the odd channels (5, T) = (0, 0) and (1, 1) 
and centrifugal potential for / = 1 (thin line). R.h.s.: Corresponding cor- 
relation functions R'^^r), R^{r), R°^'^{r) in terms of R+{r) — r. 

spin-isosp in-independent case there are now four correlation functions R^{r) 
which have to be determined. For the ATS3M potential [18,19] there are only 
three because the odd channels {S, T) = (0, 0) and (1, 1) have the same purely 
repulsive radial dependence. The {S,T) = (0,1) and (1,0) potentials have a 
repulsive core and an attractive tail so that the procedure discussed in 2.4 is 
used to construct a correlation function of finite range from eq. (82) which is 
then parameterized in the form (96). The resulting parameters are given in 
table 3 together with the correlation volumes. The correlation functions are 
shown on the r.h.s. of fig. 12. 

For the repulsive channels {S, T) = (0, 0) and (1, 1) with an odd parity in 
the relative wave function the exact correlation function obtained by eq. (82) 
does not provide a unique scale parameter j3. As can be seen from fig. 12 
the centrifugal barrier for I = 1 states dominates the repulsive potential for 



channel 


a [fm] 


/3 [fm] 


V 


Vc [fm^] 


even 


S = 0,T = 


1 


1.81 


1.07 


0.67 


0.43 


even 


S = 1,T = 





1.43 


0.95 


0.78 


0.21 


odd, 


S = 0, T = 





2.30 


1.00 


0.90 


0.004 


odd, 


S = l, T = 


1 











Table 3 

Parameters of the different correlation functions R^{r) and corre- 
sponding correlation volumes Vc- 
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[AMeV] He 



16, 



20 



19.63 



22.37 



o 



23.81 



40, 



Ca 



27.72 



14.89 



1.38 



17.68 



2.53 



18.48 



1.82 



21.21 



1.74 



48Ca 
27.05 



20.99 



2.31 



-6.97 



<-7.2 -5.11 



-13.51 



-20 



-15.15 



-16.67 



-6.74*''"^'^^ -8^<--8.65 -7.04 



<-7.57 



-26.61 



-27.48 



-30.56 



-36.09 



-34.09 



Fig. 13. Uncorrelated and correlated kinetic, potential and total energies 
for the modified Afnan-Tang S3 potential, otherwise same as fig. 8. 
Arrows indicate the results of a Yakubovski calculation for ^He [20], of 
a FHNC calculation for ^^O [21] and of CBF calculations for 40,48ca 
[22]. 



distance r > 0.7fm and hence even the uncorrelated / = 1 state u{r) = 
rji{\/mE r) / -y/iVo r^/ \/Nq has already a small amplitude in the repulsive 
region. Therefore, a rather small correlation function is needed to correct the 
small distance behaviour of the repulsive wave function, rj which controls the 
correlation function at small r is chosen to minimize the repulsion in the two- 
body system, while the i^O ground state was employed to choose a and f3. 
The resulting parameters are given in table 3 and the corresponding R°^'^{r) = 
R^{r) = R]^{r) is depicted on the r.h.s. of fig. 12. 



4-2 y4 = 2 ■ ■ - 48 nuclei with the modified Afnan-Tang force 



The correlation functions given by the parameters of table 3 are now used to 
calculate the energies of several nuclei. The results for "^He, ^^C, ^^O, ^''Ca and 
"^^Ca are summarized in fig. 13 in terms of the uncorrelated and correlated 
kinetic and potential energies as already explained in connection with fig. 8. 
The total energy ( H'-^^ ) is compared to results of different methods such as 
Yakubovsky calculations [20], FHNC [21] and CBF [22]. 

Considering the simplicity of the uncorrelated state and the fact that the cor- 
related two-body Hamiltonian is the same for mass number A = 2 up to 
A = 48 the agreement with much more elaborate methods is striking. For de- 
formed nuclei like i^C, ^°Ne or ^'^Si the variational many-body calculations are 
very difficult in terms of computing. With the new Unitary Correlation Op- 
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X-'-int 1 


( V ) 


{H ) 


exact energy 


/ iT'^^ \ 1/2 

{Rs ) ' 




12.33 


-14.50 


-2.17 


-2.22 


1.93 


^He 


78.5 


-106.4 


-27.9 


-28.8 [20]; -24.2 [19] 


1.35 




99.3 


-113.7 


-14.4 




1.97 


12 Q 


268.5 


-329.8 


-61.3 




2.36 


160 


381.0 


-488.9 


-107.9 


-105.3 [21]; -107.7 [19] 


2.28 


19p 


448.2 


-558.8 


-110.6 




2.52 




501.1 


-628.4 


-127.4 




2.56 


24Mg 


606.7 


-757.6 


-150.9 




2.74 




678.2 


-851.1 


-172.9 




2.82 




705.4 


-896.1 


-190.7 




2.88 




1108.7 


-1443.6 


-334.9 


-346.0 [22]; -335.6 [19] 


2.93 


4«Ca 


1298.4 


-1636.1 


-337.7 


-363.36 [22] 


3.20 



Table 4 

Results of a variational calculation with a single Slater determinant 
(except for ^H) and a fixed correlator for several nuclei with the ATS3M 
potential. Energies are given in MeV, radii in fm. 



erator Method combined with the uncorrelated states of Fermionic Molecular 
Dynamics [23] deformed nuclei do not pose special problems. For example for 
^°Ne the minimization takes about 3 minutes CPU time on an IBM RS6000, 
for '^'^Ca about one hour. 

In table 4 ground state energies of a variety of nuclei, including those of fig. 13, 
are listed and compared to literature values if available. 

Considering the fact that for all nuclei the energy results from the subtraction 
of a large positive kinetic energy and a similarly large potential energy, which 
are about 4 respectively 5 times larger than the total energy, it is surprising 
that the outcome of the present method lies within the uncertainties of "ex- 
act" many-body calculations. One reason is probably that the probability in 
the two-body density which is moved out of the repulsive core region is put 
where the potential is most attractive, so that a large part of the attractive 
correlations are also taken care of. Otherwise it is hard to imagine that a single 
Slater determinant which represents only the mean-field behavior is sufficient 
as an uncorrelated state. 
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The root-mean-square radii of the matter distribution which arc given in the 
last column of table 4 are calculated with the correlated operator in two-body 
approximation. 



Rf, 



r2 

"cm 



(109) 



-C2 



i i<j 



where the centre of mass position operator 



A 



E^(0 



(110) 



commutes with the correlator C. The expression (109) is a straightforward 
result of the unitary correlator method explained in section 2. 

The contribution {Rg ) of the second order to ( i?s ) is in all cases only 
in the last digit which is displayed. This shows that the matter distribution in 
coordinate space is changed very little by the repulsive two-body correlations. 



5 Summary 



A new concept, the Unitary Correlation Operator Method (UCOM), has been 
developed to describe short range correlations brought about by the repulsive 
part of the interaction. 

The unitary correlator reduces the probability to find two particles in the 
classically forbidden region of the repulsive core by shifting them away from 
each other. The relation between uncorrelated and new correlated distances 
are expressed in terms of the coordinate transformation r — > R+{r). With 
this key quantity, which may depend on spin and isospin, analytic expressions 
for the correlated states and the correlated kinetic and potential energies are 
derived. Due to unitarity correlated operators can be evaluated easily, for 
example a local potential v{r) between two particles transforms as v{r) c ~ 
v{r) c — v{c~^ rc)= v{R+{r)). 

For a many-body system the physical challenge is to find an optimal -R+(r) 
which describes well the depletion in the two-body density at short distances, 
which cannot be represented by the uncorrelated states. Due to unitarity this 
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reduction in probability results in an increase at larger distances. If the two- 
body potential has an attractive area outside the repulsion it is natural to 
enhance the probability there. These type of repulsive-attractive potentials 
provide an inherent length scale for the range of the correlator. 

For numerical ease one wants to use the two-body approximation for corre- 
lated operators and states and hence the correlation range should be smaller 
than the mean particle distance so that the probability to find two particles 
within the correlation volume Vc is small and the chance that three particles 
are simultaneously inside Vc is negligible. If this condition is not fulfilled be- 
cause the density is too high, the correlation range in which i?+(r) differs from 
r has to be reduced and at the same time the flexibility of the uncorrelated 
states has to be increased to include medium range correlations. This is an 
important aspect of the method: the unitary correlator and the uncorrelated 
states have to harmonize. 

For nuclei (at least without tensor forces which we have not investigated yet) 
this seems to be the case because the calculations presented show that a 
single Slater determinant plus a suitably chosen unitary correlator reproduce 
the "exact" ground-state energies of nuclei between A — 2 and 48 within a 
few percent. 
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A Appendix 

A.l Malfliet-Tjon V potential 
The MTV potential [13] is given by 



^ exp{-/3ir} 
1=1 ' 



i 


7i [MeV fm] 


Pi [fm-i] 


with I 


-584 


1.55 


2 


1458 


3.11 



A. 2 Modified Afnan-Tang S3 potential 

The modified Afnan-Tang S3 potential [19] is given by 

V = Vo{r) IIoi + vi{r) Uiq + Vodd{r) {Uqq + rtn) with 

3 

'^s(r) = Yl 75,j exp {-(3s,i r^} 
1=1 

and the parameters 



(A.l) 



(A.2) 





S 


= 


S 


= 1 


odd 


i 


ls,i [MeV] 


(is,^ [fm-2] 


7S,, [MeV] 


Ps,i [fm-2] 


75,i [MeV] ps,i [fm-2] 


1 


1000 


3.0 


1000 


3.0 


1000 3.0 


2 


-166 


0.8 


-326.7 


1.05 




3 


-23 


0.4 


1000 


0.6 





Table A.l 

Parameters of the modified Afnan-Tang S3 potential. 

The operators II^t = 11^ project on two-body states with total spin S 
and total isospin T, respectively. 
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